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ABSTRACT 

Secondary anisotropics in the cosmic microwave background are a treasure-trove of cosmological information. 
Interpreting current experiments probing them are limited by theoretical uncertainties rather than by measurement 
errors. Here we focus on the secondary anisotropics resulting from the thermal Sunyaev-Zel'dovich (tSZ) effect; 
the amplitude of which depends critically on the average thermal pressure profile of galaxy groups and clusters. To 
this end, we use a suite of hydrodynamical TreePM-SPH simulations that include radiative cooling, star formation, 
supernova feedback, and energetic feedback from active galactic nuclei (AGN). We examine in detail how the 
pressure profile depends on cluster radius, mass, and redshift and provide an empirical fitting function. We employ 
three different approaches for calculating the tSZ power spectrum: an analytical approach that uses our pressure 
profile fit, a semi-analytical method of pasting our pressure fit onto simulated clusters, and a direct numerical 
integration of our simulated volumes. We demonstrate that the detailed structure of the intracluster medium and 
cosmic web affect the tSZ power spectrum. In particular, the substructure and asphericity of clusters increase 
the tSZ power spectrum by 10-20% at £ ~ 2000-8000, with most of the additional power being contributed by 
substructures. The contributions to the power spectrum from radii larger than Z?soo is ~ 20% at £ = 3000, thus 
clusters interiors (r < R500) dominate the power spectrum amplitude at these angular scales. 
Subject headings: Cosmic Microwave Background — Cosmology: Theory — Galaxies: Clusters: General — 
Large-Scale Structure of Universe — Methods: Numerical 



1. INTRODUCTION 

As cosmic microwave background (CMB) photons travel 
through the diffuse hot gas comprising the bulk of baryons in 
C"} I galaxy clusters, a fraction of them are upscattered by the gas 
Q\ • in a process called the thermal S unyaev-Zerdovich (tSZ) ef- 
f^i ' feet (ISunvaev & Zeldov ich 1970). This scattering produces a 
y—i ' unique spectral signature in the CMB, with a decrement in ther- 
7— I . modynamic temperature below v ~ 220 GHz, and an excess 
k/ ' above. The tSZ effect is typically seen on arc-minute scales, 
. £^ \ and is referred to as a secondary anisotropy, as it originates 
, between us and the surface of last scattering, unlike the pri- 
j_j ■ mary CMB anisotropies. In the non-relativistic limit, the tSZ is 
d ' directly proportional to the integrated electron pressure along 
the line-of-sight. It typically traces out the spatial distribu- 
tion of clusters and groups, since the hot intracluster medium 
(ICM) dominates the line-of-sight pressure integral. Thus, the 
tSZ provides an excellent tool to examine the bulk of cluster 
bary ons. Found at the intersections of filaments in the cosmic 
web dBond et al.l 19961) . clusters form at sites of constructive in- 
terference of long waves i n the primordial density fluctuations, 
the co herent peak-patches (Barde en et al.ll9 86: Bond & Myers 
119961) . Clusters are sign posts for the growth of structure in the 
Universe, and are a potentially powerful tool for probing un- 
derlying cosmological parameters, such as w, the dark energy 
pressure-to-density ratio. 

The angular power spectrum of the tSZ effect is extremely 
sensitive to cosmological parameters like erg, the root mean 
square (RMS) amplitude of the (linearized) density fluctuations 
on 8/T 1 Mpc scales. In fact, the amplitude of the tSZ power 
spectrum scales at least as steeply as the seventh power of erg 



dBond et all 120021: iKomatsu & Seliakll200l iBond et al.l 12005: 
iTrac et al.l201 ll) and improving the constraints on erg will aid in 
breaking the degeneracies found between erg and w when using 
only primary CMB constraints. An advantage of using the tSZ 
angular power spectrum over counting clusters is that no ex- 
plicit measurement of cluster masses is required. Also, lower 
mass, and therefore fainter, clusters that may not be signifi- 
cantly detected as individual objects in CMB maps contribute 
to this statistical signal. However, disadvantages of using the 
tSZ angular power spectrum include potential contamination 
from point sources and that no redshift information from the 
clusters is used. 

Previous observatio ns by the Be r keley- Illinois-Maryland 
Association (BIMA, Daw son etaTI 120061). the Atacam a 
Path-finding Experiment (APE X-SZ. iReichardt etai]|2009bl) . 
the Quest at DASI (QUaD, iFriedman et al.l 120091) . Arc- 
minute Cosmology B olometer Array Receiver (ACBAR, 
Reichard t et al.l l2009al). a nd the Cosmic Background Imager 
(CBI, ISievers et al.l 120091) all measured excess power above 
that expected from primary anisotropies, which have been 
attributed to some combination of the tSZ effect and point 
source contamination. The measurements from these experi- 
ments provided upper limits to the tSZ power spectrum am- 
plitud e. More recently, the Atacama Cosmo logy Telescope 
(ACT, iFowler et al.ll2010t iDunkley et al.ll2010l) and the South 
Pole Telescope (SP T. lLueker et aLlboiOtlshirokoff et al1 l20 10; 
iKeisleret all 1201 1|) have detected the SZ effect in the CMB 
power spectrurrfl The results from ACT and SPT emphasize 

1 The Planck collaboration has released some early SZ science (e.g., 
Planck Collaboration et al. 2011a b c), but to-date there have been no power 
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that the "sweet spot" for measurin g the tSZ signal is between 
I ~ 2000-4000. Silk damping (ISilkl 119681) suppresses the 
power of primary anisotropies so that their contributions to the 
power spectrum are much smaller than the tSZ contribution at 
even higher I. At these scales there are important additional 
contributions to the power spectrum from the kinetic SZ (kSZ) 
effect, which arises from motions of ionized gas with respect 
to the CMB rest frame, as well as dusty star-forming galaxies 
and the radio galaxies, both of which appear as point sources. 
All these signals increase the uncertainty when determining the 
tSZ power spectrum, and hence the parameters derived there- 
from. 

Three main tools have been used to estimate the tSZ power 
spectrum: Analytic models, semi-analytical models, and 
numerical simulations. They have been used to derive several 
different templ ates for the predicted tSZ power s pectrum (e.g., 
Cole &Kaiserl [l988t IMakino & Sutol Tl99l (da Silva et al l 
20001: iRefreeier et alJ l2000t iHolder & Carlstroml 120011: 
Zhang & Pen! 120011: ISpringel et all 1200 ll iKomatsu & Seliak 
200 21; IZhangetalJ 120021: iBondet all 120051: ISchafer et al. 
200 6al["bl; iBattaglia et all l2010t IShaw et alJ 120101: iTrac et al. 
2011; Efstathiou & Migliaccio 2011). There are both shape 
and amplitude differences between these three approaches that 
compute the tSZ power spectru m; comparisons are required 
to understand these differences (Refregieret al. 2000). At the 
base of these differences is the cluster electron pressure profile, 
since it is a crucial and uncertain component in the analytical 
thermal SZ power spectrum calculation. The electron pressure 
profile is directly related to the total thermal energy in a cluster 
and is sensitive to all the complicated gastrophysics of the 
ICM. For example, some of the ICM processes that should 
be included are radiative cooling, star-formation, energetic 
feedback from AGN and massive stars, non-thermal pressure 
support, magnetic fields, and cosmic rays. Deviations from 
an average pressure profile, i.e., cluster substructure and 
asphericity will also contribute to the tSZ power spectrum. But 
how much? 

The inclusion of AGN feedback is vi tal to any tSZ power 
spectrum template (Batta glia et al.ll20lol) . Furthermore, an en- 
ergetic feedback source (AGN feedback being the most pop- 
ular) seems to be an important addition to any hydrodynam- 
ical simulation, since simulations with only radiative cooling 
and supernova feedback have problems with exc essive over- 
cooling in cluster centers (e.g. jLewis et alj|2000i) . This over- 
cooling results in too many stars being produced out of ICM 
gas reservoir, which alters the thermal and hydrodynamic struc- 
ture of ICM in a way that is inconsistent with observational 
data. 

In this paper we present a detailed comparison of the three 
approaches used to calculated the thermal SZ angular power 
spectrum. This comparison allows us to identify and quantify 
the differences between each method. Section|2]briefly summa- 
rizes the simulations used in this work and Section [5] outlines 
the calculation of the analytical tSZ angular power spectrum. 
In Sections H] and we present our results for numerical aver- 
age thermal pressure profiles and a detailed analysis of the tSZ 
power spectrum, respectively. In Section|6]we provide updated 
constraints on a% using the new ACT and SPT measurements 
of the CMB power spectrum at high £, and we summarize our 
results and conclude in Section [7] 



spectrum results. 



2. COSMOLOGICAL SIMULATIONS AND CLUSTER DATA SET 

We use a modified version of the smoothed parti cle hydro- 
dynamical (SPH) code GADGET-2 (Springel 2005) to simu- 
late cosmological volumes. We use a suite of 10 simulations 
with periodic boundary conditions, box size 165/r'Mpc, and 
with equal numbers of dark matter and gas particles A^dm = 
Ngas = 256 3 . We adopt a flat tilted ACDM cosmology, with total 
matter density (in units of the critical) fl m = f^DM + ^b = 0.25, 
baryon density = 0.043, cosmological constant S1a = 0.75, 
a present day Hubble constant of Hq = 100/? km s~ Mpc~ , 
a scalar spectral index of the primordial power-spectrum n s 
= 0.96 and <j% = 0.8. The particle masses are then m gas = 
3.2 x IO^-'Mq and m DM = 1.54 x IO^/z^Mq. The mini- 
mum gravitational smoothing length is e s = 20 h kpc; our SPH 
densities are computed with 32 neighbours. 

We include sub -grid models for AGN feedback 
(Battaglia et alJ l2010h . r adiative cooling , star formation, 
and supernova feedback dKatz et alJ [l996t lHaardt & M adau 
Il996t ISpringel & Hernquistl I2003IT The AGN feedback 
prescription included in the simulations (for more details see 
Battaglia et al] 120101 1201 ll) allows for lower resolution and 
hence can be applied to large-scale structure simulations. It 
couples the black hole accretion rate to the glo bal star forma- 
tion r ate (SFR) of the cluster, as suggested by Thom pson et alJ 
(2005). The thermal energy is injected into the ICM such 
that it is proportional to the star-formation within a given 
spherical region. Throughout this work we will refer to these 
simulations as AGN feedback. 

We adopt the standard working definition of cluster radii Ra 
as the radius at which the mean interior density equals A times 
the critical density, p a (z) (e.g., for A = 200 or 500). For clarity 
the critical density is 

Pcr(z) = ^[n m (i+z) 3 +nA\. (1) 

Here we have assumed a flat universe + Cl\ = 1) and are 
only interested at times after the matter-radiation equality, i.e., 
the radiation term with il r is negligible. It is important to note 
that all masses and distances quoted in this work are given rel- 
ative to h = 0.7, since most observations are reported with this 
value of h. 

3. THE ANALYTIC CALCULATIONS OF TSZ ANGULAR POWER 

SPECTRUM 

The tSZ can be adequatel y modelled as a rando m distributed 
Poisson process on the sky (Cole & Kaiser 19880- There are 
two components in this model that are required for a statistical 
representation of the secondary anisotropies: (1) The number 
density for objects of a given class; and (2) the profile of the 
same object and class, centered on its position. We focus on 
groups and clusters, since they are the dominant source of tSZ 
aniso tropies. This approac h is referred to as the halo formalism 
(e.g. JCole & Kaiserlll988l) . 

The non-relativistic tSZ signal is the line-of-sight integration 
of the electron pressure, 

^ = my = m— 2 1 p e (w , (2) 

1 m e c J 

2 Note that we are not including the contributions from the clustering of 
cluste rs, since this is sub-dominant on scales of I > 300 I Komatsu & Kitayama 
[1991 . 
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where f(v) is the spe ctral function for the tSZ 
(ISunvaev & Zeldovichl 119701) . y is the Compton-y param- 
eter, or is the Thompson cross-section, ra e is the electron mass 
and P e is electron pressurfl For a fully ionized medium, the 
thermal pressure P t h = P e (5X H + 3)/2(X H + 1) = 1.932P e , where 
X H = 0.76 is the primordial hydrogen mass fraction, and P t j, is 
the thermal pressure. 

We adopt the successful analytical ansatz for halo number 
density as a function of mass 



dn{M,z) 



R(M) da(M,zf 



dM 2M 2 3a(M,z) 2 dR(M) 



f(a(M,z)) (3) 



where a(M,z) is the RMS variance of the linear density field 
smoothed on the scale of R(M), and f(a) is a functional form 
determined from N-body simulations (e.g., Henkins et aDl2001b 
IWarren et al.ll2006T:lTinker et alj|2 008). In this work we use the 
mass function from lTinker et al.l (12008b for the analytic calcu- 
lations. Note that the tSZ power spectrum is only mildly sensi- 
tive to the particulars of the mass function (Komatsu & Seliak 
120021) . 

The tSZ angular power spectrum at a multipole moment £ is 



Q.tsz = fiyY 



(4) 



where y^(M,z) is the form factor, which is proportional to the 
Fourier transform of the projected electron pressure profile, P e . 
We do not include h igher order relativistic corrections to f(u) 
dNozawa et al.ll2006h . 

The functional form of yi(M,z) can be determined em- 
pirically in observations or simulations (e.g., iNagai et aD 
120071: lArnaud etal.1l 2010). or can be determined an alytically 
(e. g.. komatsu & S eliak 2001; Ostriker et ail 120051) . Follow- 
ing Komatsu & Seliak (2002) we compute y~i(M,z) assuming 
spherical symmetry and using Limber's approximation, 



n(M,z) = 



Airr s <Jj 



m e c- 



2d , .sin(£x/4) 

x w ^ilT^ 



(5) 



where x = r/r s is a dimensionless radius, £ s = D&/r s is the 
corresponding angular wave number, and Da is th e angular di- 
ameter distance. We follow Navar roet al.l ( 1 19971) in our defi- 
nition of the scale radius in a cluster with concentration cnfw, 
= ^vir/cNFw. H ere we use a fitting formula for cnfw from 
Duf fy et all (120081) and t he definition for the virial radius from 
Bryan & NormarT(fT998l) . 



3M V 



4 7 rA cr (z)p cr (z) 



1/3 



(6) 



where A cr (z) = 187r 2 + 82[rj(z)- l]-39[f2(z)- l] 2 and Q(z) = 

n m (i+z) 3 [n a (i+z?+n A y\ 

The dominant source of uncertainty in Q. t sz comes from 
y~i(M,z), since one can easily calculate the volume element for 
a given cosmology , and the mass function is known to 5 - 10% 
(Tinkeretalj2008). Thus, the pressure profile is the critical in- 
put into the analytical tSZ angular power spectrum. We would 
ideally like to know y~e(M, z) as well as we know the mass func- 
tion. This requires an understanding of the detailed physical 
processes which affect cluster pressure profiles. 

3 Here we have ignored the temperature of the CMB, 7cmb. since Tqmb *C 
T e , hence n e &B (7s - Tqmb) — «c k% T e = P e . 



The Gaussian and non-Gaussian variances of the power spec- 
trum are also calculated using the halo formalism (Bond 1996; 
| Cooravll200ll; iKomatsu & Seliak! Hool IZhang & Shethl 120071: 
Sha w etal.1 120091) . again neglecting the clustering of clusters 
term. The full-sky variance is 



'.tSZ ■ 



2(Q.tsz) 2 



2l+\ 



4tt 



(7) 



where Tu> is the trispectrum; see Equation The variance is 
inversely proportional to the sky area covered, so for a fraction 
/sky of the sky covered, a 2 u , tSZ oc 1/ / s k y . In this work we will 
present the diagonal part of the covariance; the diagonal of the 
trispectrum is 



T m sz = f(v) 4 J ^dz J 



^lUM^dM. 



(8) 



4. THE THERMAL PRESSURE PROFILE 



The cluster thermal pressure profile is the most uncertain 
component of the tSZ power spectrum. In this section we 
use a large sample of clusters from hydrodynamical simula- 
tions and explore the mean cluster pr ofile and the subtle differ - 
ences from self-similar scaling (e.g.. lKaiserlll986fc IVoitll2005b . 
Comparisons between the latest pressure profiles from analyt- 
ics, observations, and simulations have sh own that they are in 
reasonable agreement with one another (lArnaud et all [2010; 
IShaw et alJl20Tot iTrac et al.ll20TTt ISun et al.ll2011l) . Previous 
work has shown that AGN feedback can alter the pressure pro- 
files, though the profil es are comparable t o previous simula- 
tions and observations (Battaglia et al. 2010). We show the de- 
pendence of the pressure profile on the cluster mass and red- 
shift and explore deviations from the self-similar scaling. 

4.1. Fitting Pressure Profiles from the Simulations 

We apply the following four-step algorithm to compute the 
average thermal pressure profiles in our simulations. First, 
we find all clusters in a given snapshot using a friends-of- 
friends (FOF) algorithm dHuchra & Gellerll 19821) using a link- 
ing length of 0.2 and an M FO f mass cut of 1.4 x 10 I3 M Q . 
Second, starting with a position and radius derived from the 
FOF results, we find the final cluster positions by recursively 
shrinking the radius of the sphere examined, and re-center on 
its center of mass. Given the cluster center, we then calcu- 
late the spherical-overdensity mass and radius, Ma and Ra- 
Third, we calculate the thermal pressure profile for the en- 
tire sample of clusters in spherical shells, with the shells de- 
fined relative to R& (for the pressure profiles, we use A = 200). 
To facilitate profile comparisons and cluster stacking, we nor- 
malize each profile by the sel f-similar amp l itude for pr essure 
P A = GM A Ap cr (z)/b/(2P A ) (lKaiserlll986t fVoit 2005]), with 
/b = fib/ S^m- Finally, we form a weighted average of these pro- 
files by stacking clusters in a given redshift and mass bins. We 
use the integrated Compton-y parameter as our weighting func- 
tion, 



P e (r)47rr 2 dr ocE th «/? A ), 



(9) 



Jo 

The stacked average profiles P^ = (P±/Pa) are then fit to a 
restricted version of the generalized NFW profile, 



Pbt = Po {x/x c f [1 + (x/x c ) a ]~ P , x = r/R A: (10) 
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Figure 1. The normalized average pressure profiles and parametrized fits to these profiles from simulations with AGN feedback scaled by (r/i?2oo) 3 . m mass bins 
(left panel) and redshift bins (right panel). Here we have independently fit each mass and redshift bin. The grey band shows the standard deviation of the average 
cluster in the most massive bin (left) and lowest redshift bin (right). In both panels we illustrate the radii that contribute 68% and 95% of the total thermal energy, 
7, centered on the median, by horizontal purple and pink error bars. The bottom panels show the percent difference between the fits and the average profiles. The 
generalized NFW profile with fixed a and 7 fits the average profiles well in the majority of the mass and redshift bins, with deviations within ~ 5% of the mean. The 
upturns at large radii are due to contributions from nearby clusters and substructure. 




Figure 2. The normalized average pressure profiles and constrained fits to these profiles from simulations with AGN feedback scaled by (r/i?2oo) 3 , f° r mass bins 
(left panel) and redshift bins (right panel). The constrained fit is a global pressure profile, as described in Section l4^2l with parameters in Table[T] It differs from the 
fits in Figure [T] where each bin was fit independently. The grey band shows the standard deviation of the average cluster in the most massive bin (left) and lowest 
redshift bin (right). In both panels we illustrate the radii that contribute 68% and 95% of the total thermal energy, Y, centered on the median, by horizontal purple and 
pink error bars. The bottom panels show the percent differences between the constrained global fits and the average profiles. The constrained fits match the average 
profiles well in the majority of the mass and redshift bins and the deviations are within ~ 10% of the mean. The upturns at large radii are due to contributions from 
substructure and nearby clusters. 
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where the fit parameters are a core-scale x c , an amplitude Po 
and a power law index j3 for the asymptotic fall off of the pro- 
file. There is substantial degeneracy between fit parameters, 
so we fix a = 1 .0 and 7 = -0.3 (as suggested by Nagai et a"T] 
12007b lArnaud et alJl2010h . We find that fitting for all parame- 
ters did not provide a significantly better fit than when a and 
7 were fixed. However, without fixing a and 7, a direct com- 
parison of fit parameters between different mass and redshift 
slices was not meaningful. We find the best-fit parameters us- 
ing a non-linear least squares Lev enberg-Marquardt approach 
(iLevenbergl 1 944t iMarquardt l 1 9631) . We weight each radial bin 
by the internal variance of the cluster profiles within that bin. 

In FigureQ] we show the mass and redshift dependence of the 
average cluster thermal pressure profile and the corresponding 
parametrized fits to these profiles. We scale the pressure pro- 
files by jc 3 , such that the height corresponds to the contribution 
per logarithmic radial interval to the total thermal energy con- 
tent of the cluster (cf. horizontal purple and pink error bars 
for the radii that contribute 68% and 95% of the cluster ther- 
mal energy). In the bottom panels of Figure Q] we highlight 
the residuals from the smoothed fitting function by showing 
the relative difference in per cent, Af t h =100 (P&t~ Ah) /P±- 
The fitting function, Equation [10] provides an accurate fit over 
all mass and redshift ranges, with a majority of the deviations 
from the average profile being < 5%. 

We find that there are subtle dependencies on the cluster 
mass and redshift (cf. Table [U, which suggests that neglecting 
these dependencies would not yield the required 5- 10% pre- 
cision needed for calculations of tSZ power spectrum. We also 
find that there are contributions to the average pressure pro- 
file at larger radii from substructure and nearby clusters, which 
cause relative deviations from the mean profile > 5%. In a 
companion paper, we also show that substructure affects the 
kinetic support in clust er outskirts and the shape of the ICM 
shape at similar radii (Battaglia et al.ll201 lb . In these regions 
(redshift dependent, but typically > 2/?20o) P&t often deviates 
from P± by more than 5%. We chose not to model this behavior 
because of two reasons. First, the problem of double-counting 
SZ flux: the large volume contained within the radius that con- 
tains 95% of the total SZ flux, r < 4R 200, necessarily leads 
to overlapping volumes of neighboring clusters, especially at 
high-redshift. Second, the total SZ flux of an increasing pres- 
sure profile, scaled by x 3 , does not converge and an arbitrarily 
chosen radial cutoff would substantially contribute to the re- 
sulting power of the tSZ power spectrum. Because we weight 
by the variance within radial bins, these contaminated regions 
are naturally down- weighted in the profile fits. 

4.2. Constrained Thermal Pressure Profile Fits 

In this section we derive a global fit to our pressure profiles 
as a function of mass and redshift. We find treating each pa- 
rameter as a separable function of mass and redshift gives good 
results, with the fit parameters constrained to be of the follow- 
ing form: For generic parameter A, we have 

For each of Po, j3, and x c , we find a m by fitting to the z = snap- 
shot, and we find a z by fitting to clusters with 1 . 1 x 1O 14 M < 
M20Q < 1.7 x 10 14 M Q . The weights used in the fits were the 
inverse variance of the fit parameters when fitting each indi- 
vidual cluster in that mass/redshift bin. With these fit parame- 
ters, presented in Table Q] and using Equations ( fTOb and ( fTTI ). 



Table 1 

Mass and Redshift Fit Parameters from Eqns. 4 1 0t and (TTJ. 



Parameter 


Am — A z 


Otm 




Po 


18.1 


0.154 


-0.758 


Xc 


0.497 


-0.00865 


0.731 




4.35 


0.0393 


0.415 



The input weights are chosen to be the inverse variances of 
fit parameter values from the individual pressure fits for each 
cluster within the bin. 

we now have a global model for the average electron pressure 
as a function of cluster radius, redshift, and mass. Hereafter 
we refer to this global empirical description as the constrained 
pressure profile. In Figure [2] we compare the constrained fits 
to the stacked averages. With fewer degrees of freedom, the 
constrained fits will naturally not be as accurate as fitting each 
mass/redshift bin completely independently, but we find that 
the mean recovered profile is accurate to 10% and corresponds 
well to the accuracy with which we intend to measure the tSZ 
power spectrum. 

The average of t his global constraine d pressure profile at 
z = (as reported in Battagli a et alj|2010b compares well with 
the average universal pressure profile from a representative 
XXM-Newton sample of nearby s ystems for th e region r < Pgoo 
where X-ray data is ava ilable (Arnau d et al.ll2010b . We defer 
the reader to Battagl ia et al.l ((2010) for a detailed discussion 
and comparison to other numerical and observational work. 
However, we stress that the global constrained pressure pro- 
file of Equation (fTTT i models the mass dependence and predicts 
a redshift evolution that shows small but noticeable deviations 
from the self-similar scaling on account of the radiative gas 
physics including AGN feedback. This has clear implications 
when analyzing SZ measurements of non-local clusters. 

In Figure [3] we present projected 30 GHz temperature maps 
of 4 sample clusters (cut at a spherical radius of 6P500), their 
expected maps from the global constrained fit, and the errors 
in the predicted temperature. A quantitative co mpa rison of the 
tSZ power spectrum is deferred until Section 15.11 Hereafter, 
we refer to the predicted temperature maps as pasted profile 
maps. Note that this is not a representative sample of the dif- 
ference between the pasted profiles and the simulations. In- 
stead, we attempt to show different size clusters across differ- 
ent redshifts and illustrate the scales of the deviations from the 
constrained fit, primarily resulting from substructure and mis- 
centering, since the cluster center of mass does not necessarily 
line up with the peak of the projected pressure. In the rightmost 
panel of Figure [3] we show the residuals amplitudes between 
the simulated cluster projections and the pasted profile from the 
constrained fits. We find that these profiles are within ~ 10% of 
the actual simulated cluster, which is similar to the differences 
found in the bottom panels of Figure[2] These substructures are 
significant on scales of tens of arc minutes for nearby massive 
clusters and scales of arc minutes for higher redshift clusters, 
corresponding to £ ~ 1 000 - 1 0000. 

4.3. Analytic Assumptions in the Thermal Pressure Profile 

Analytic and semi-analytic models typically rely on assum- 
ing an pressure-density (P—p) relation and some form of hy- 
drostatic equilibrium (HSE), possibly incl uding non-thermal 
support terms. Fully an alytic models, (e.g., Komatsu & Seljak 
20021 IShaw et alj |2010). apply HSE to theoretical, spherically 
symmetric dark matter potentials. Semi-analytic models, (e.g., 
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Figure 3. A comparison of four projected pressure maps of simulated clusters to the projected pasted-profile maps. From left to right, the panels show the simulated 
clusters (cut at a spherical radius of 6Rjoo)> the projected pasted profiles from the constrained fit, and the difference map between the two. The maps show the 
temperature decrement —AT in units of fiK, at a frequency of 30 GHz. The difference maps, ST, illustrate the scales and amplitudes of the residuals between the 
simulated clusters and the projected pasted profiles. Note the color scale is logarithmic for the left two panels (from -0.1 fiK to -300 fiK), while it is linear for the 
difference map (from — 30/jK to 30/^K). For all panels the left and top axes are in units of Mpc and the bottom and right axes are in units of arc-minutes. 
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Figure 4. The assumption of a constant thermal or total logarithmic slope of 
the P-p relation, T,), or r to t, as most analytic models assume, is not consistent 
with the results from our simulations. r t h (solid line) and r l0 t (dashed line) are 
shown as functions of radius at z = and z = 1 from simulations with AGN feed- 
back. For comparison, we show the total adiabatic index used by Shaw et al. 
(2010), and we find that the differences increase at larger radii, especially at 
high redshifts. 

lSehpaletaill2010t iTrac et al.ll201 lb . take dark matter simula- 
tions, and paste baryons on top of the dark matter potential 
wells, again using (possibly corrected) HSE and an P-p re- 
lation. The results from both classes of models, then, rely 
critically on the input P-p relation and are sensitive to depar- 
tures from HSE. In contrast, empirical fits to the average cluster 
pressure profile derived from simulations have a key advantage 
over analytical models because the simulations naturally ac- 
count for kinetic pressure support from non-thermalized bulk 
flows which provide substantial support in the outer parts of 
clusters, but do not contribute to the tSZ. They also make no as- 
sumptions about HSE (which is grossly violated during, for in- 
stance, mergers), and rather than forcing an P-p relation, they 
track the flow of energy into and out of the ICM. 

The (semi-)analytic calculations cast the P-p relation in 
terms of a pressure law P oc p 1 ', and usually assume a constant 
F, where P can be either the thermal pressure P t h which is the 
source for the tSZ effect, or the total pressure, P m = P t h + Pnt, 
where P nt is any non-thermal support, principally kinetic mo- 
tion of the IQVQ The total pressure is the input to the equation 
of hydrostatic equilibrium that reads for spherical symmetry 

dP tot /dr = -GM«r)p/r 2 . (12) 

We present the effective adiabatic index, V = dlogP/dlogp, 
as a function of cluster radius in Figure [4] We find that the as- 
sumption of constant T is grossly violated, particularly in the 
outer parts of clusters, and for P t h. These results stress the im- 
portance of deriving pressure profiles from observations and 
hydrodynamical simulations, particularly as good-quality ob- 
servational data from cluster outskirts is in short supply. 

5. THE TSZ POWER SPECTRUM IN DETAIL 

In this section we compare three different ways of calcu- 
lating the tSZ power spectrum: directly projecting the elec- 
tron pressure in the simulations, taking the simulation cluster 

4 Some older models have ignored kinetic support entirely, in which case 
Pm = Pth- 



catalogs and projecting our constrained global pressure profile 
onto the cluster locations (the "pasted-profile" maps), and us- 
ing a completely analytical halo calculation. For the analytic 
calculation, we use the formalism described in Section [3] and 
the constrained pressure profile from Section [4] For the sim- 
ulation and pasted-profile approaches, the thermal Compton-y 
maps are obtained by performing a line-of-sight integration of 
the electron pressure through the entire simulation box at each 
redshift output, covering z = 0.07 to z = 5. For each redshift- 
output map we compute the average power spectrum for our 
ten simulations and add these differential power spectra upB 
This procedure uses all the information within the simulation 
volume and decreases the variance of the power spectrum, es- 
pecially at low redshifts. One benefit of this technique is that 
by summing over redshift slices after taking the power spectra, 
we ignore any correlations between different redshift slices, as 
effectively happens in nature. With more traditiona l methods 
that s tack redshift slices (such as were used in Battaglia et al.1 
120101) . care must be taken that different redshift slices do not 
project the same objects to the same locations, as that induces 
artificial correlations, potentially altering the tSZ power spec- 
trum. 

In the left panel of Figure [5] we plot the tSZ power from 
our analytical halo calculation and that from the AGN simula- 
tions. For reference, we include other tSZ power spectrum tem- 
plates dKomatsu & Seliak 20021; ISehgal et al.ll2010HShaw et all 
120 10t ITrac et al.112.01 ll) . We choose the cosmological parame- 
ters for the halo calculation to match the simulations and inte- 
grate from z = 0.07 to z = 5, so that the only possible sources 
of differences are the mass function and the pressure profiles. 
There are clear differences between the analytical halo calcula- 
tion and the complete simulation maps. The main difference at 
low £s results from shot noise within the sample of simulated 
boxes, where we had more (though consistent within the ex- 
pected error) high-mass clusters than expected, but this is only 
a 6% effect in the total power spectrum (cf. Appendix). The 
differences at higher £s arise from deviations about the aver- 
age pressure profile, including effects of cluster substructure 
and asphericity. We see these variations in the residual maps 
of individual simulated cluster projections and pasted profile 
projections (cf. Fig. |3J. We further explore these differences 
in the power spectrum between the an alytic cal cula tion and the 
simulations in the following Sections 15.11 and 15.21 It is chal- 
lenging to determine the causes for all the differences between 
our c alculations and other calculations for the tSZ power spec- 
trum dKomatsu & Seli akl 120021; ISehgal et alj|2010t IShaw et all 
l20Tot ITrac et al.ll2011l) . since the thermal pressure profile we 
use is different from the ones used by the other calculations. 
However, the reasons for the differences we find between our 
three methods, will be generally applicable to the other meth- 
ods of calculating the tSZ power spectra. 

The right panel of Figure [3] show s a direct compar ison be- 
tween our analytical model and the Sha w et aTl d2010) model. 
In both calculations, the same cluster mass function was used 
and the power spectra are scaled to the same cosmological pa- 
rameters, so the differences are related to the model for the 
thermal pressure profile. We investigate the redshift integration 
limits^, but find they do not significantly affect the differences 

5 We have selected the redshifts at which we write out the simulation snap- 
shots to be the light crossing time of the simulation; hence, the total power 
spectrum is the sum of the differential power spectra. 

6 For the remainder of this paper, we use a low redshift cutoff of z = 0.07, 
so that we can directly compare our analytic calculation to the simulations. 
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Figure 5. In the left panel, we show a comparison of the current predictions for the tSZ power spectra at 150 GHz from our simulations with AGN feedback (red line) 
and the analytical calculations using the constrained pressure pr ofiles in this work (b lue line). The standard deviation among our 10 simulations is shown with a light 
grey band. We also include the semi-analytical simulations by Sehgal et al. (2010) (pink dotted line) andlT rac et al. (2011) which includes enhanced non-thermal 
pressure support (dark green dashed line) and the fully analytical calculations by Komatsu & Seliak 12002) (orange dotted line) and Shaw et al. (2010) (light green 
dashed line). The full-width half-max values appropriate for the Planck, ACT and SPT beams are also plotted. At low-<?, our two methods of calculating the tSZ 
diverge because our simulations happen to contain a large number of high mass objects driving the power up, tho ugh the exce ss is consistent with expected Poisson 
fluctuations. At high-^ the discrepancy is the result of substructure and asphericity, as demonstrated in Sections |5 1 1 and |5~2l The right panel shows a comparison 
between the current analytic calculations for the tSZ power spectra and how the power spectrum changes with the variation of the lower redshift limit of integration. 
The variance of 1 % of the full-sky power spectrum (cf . Equation 0) is illustrated by the grey bands for the highest and the lowest redshift limits of integration. 



at £ > 1000. We present the expected mean and standard de- 
viation of 1 % of a full-sky tSZ measurement as a function of 
lower redshift cutoff, and find that the low-£ variance is sub- 
stantially suppressed by raising the low-z cutoff. On the scales 
where the tSZ peaks, we find both the mean spectrum and the 
variance are only weakly affected by varying the redshift limit 
from z=0.01 to z=0.14. Similar res ults have been fou nd when 
making intensity cuts on sky maps (Shaw et al. 2009). 

We now present power spectra calculated directly from the 
simulations. In addition to projecting the full electron pressure 
from all particles, we also take advantage of the information 
from the simulation cluster catalogs. By doing this, we can em- 
ploy mass, redshift, and radius cuts to explore the dependence 
of the full tSZ power spectrum. By pasting our global pressure 
profile to locations and redshifts of simulated clusters, we can 
also explore, without having to worry about sample variance, 
the effects of using our profile instead of the full simulation 
results. 

We use the cluster catalogs described in Section |4~T1 and re- 
mind the reader that M^of is roughly equal to M200, though 
with large scatterQ Our cluster mass function becomes in- 
complete below M200 ~4x 10 13 M Q (cf. Appendix) primar- 
ily due to our Mfof cutoff in the original cluster finding of 
1.4 x 10 13 M Q , but partially due to the linking length merg- 
ing some clusters/gro ups into nearby l arger clusters at the 
10— 1 5% level (e.g Davis et al.1 119851: lB ertsching e"r& Gelbl 
[1991 iCole & Lacevlll996t iCohn & White! 120081) . For these 
reasons we examine only cluster with M500 > 4.2 x 10 13 M Q 
when we bin clusters in mass. 

In Figure [6] we show the cumulative distribution function 
(CDF) for the tSZ power for a CDF(M >,z <) at 1= 3000. The 
CDF illustrates where the relative amount of power originates 

7 For detailed work on comparing the mass definitions o f Mppp to Ma and 
the resulting halo mass catalogs from these definitions see More et al. (201 1) 
and the references therein. 
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Figure 6. Shown is the cumulative distribution function for the thermal SZ 
power spectrum as a function of mass and redshift at i = 3000. The curves 
show the lower mass and upper redshift cutoffs that contribute [25, 50, 75]% 
to the tSZ power spectrum. At I — 3000, half the power of tSZ power spec- 
trum comes from clusters with z > 0.6, and half comes from clusters with 
M 50 o < 2 X 10 14 M(7) For comparison, the dashed green lines show the semi- 
analytical results of Trac et al. (2011), which include enhanced non-thermal 
pressure support. 

at the 25%, 50% and 75% percentile levels. Half the power at 
I = 3000 comes from clusters with z > 0.6 and half originates 
from clusters with mass M500 < 2 x 10 14 Mfr> This result is in 
general agreem ent with other work dKomatsu & Seliak1 l2002; 
iTrac et al.l201 lb . We note that the particulars of these mass and 
redshift ranges are sensitive to the input modeling of the ICM. 
The comparatively low mass and high redshift of the clusters 
and groups that make up the bulk of the tSZ signal mean that 
they have not been as well studied as more massive and nearby 
objects. Thus, the tSZ angular power spectrum can provide a 
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statistical constraint on the astrophysical processes of impor- 
tance at high redshift and in low-mass clusters. 

5.1. Contribution to the tSZ Power Spectrum in Cluster Mass 

Bins 

In this subsection, we calculate the power spectrum in mass 
bins. This allows us to isolate the differences between the 
simulations, the pasted profile maps, and the analytic calcu- 
lation, as functions of cluster mass, integrating in redshift be- 
tween z = 0.07 and z = 5. We explore both, cumulative and 
differential mass bins. We consider all gas particles (or radii) 
within 6/?5oo when projecting the thermal pressure of the sim- 
ulations. Our method takes care of not double counting the 
cluster mass in overlapping volumes of close-by clusters. In 
Figure [7] we show the power spectrum broken down into cu- 
mulative (left panel) and differential (right panel) mass bins. 
The bottom panels show the relative differences, where AQ = 
100 (Q,j; m -Q,/) /Q, s / m , with Q, s ; m denoting the power spec- 
tra from the simulations and the Q., are the power spectra from 
either the projected pasted profile maps or the analytic calcula- 
tion. 

The largest deviations between our analytic/pasted profile 
spectra and the full simulations are for the highest mass 
(M500 > 4.2 x 10 14 M Q ) clusters, particularly on small angu- 
lar scales. The deviations between the pasted profiles and the 
simulations in this mass range arise from the increased level 
of substructure and asphericity in massive clusters in com- 
parison to smaller objects due to the more recent formation 
epoch of large systems in a hierarchical structure formation 
(IWechsler et a Jl2002t IZhao et aUl2009t iPfrommer etaDl2011h 
Battag lia et al 1). The high-mass difference between the 
fully analytic tSZ spectrum and the simulation results reflects 
our overabundance of high-mass clusters due to shot noise rel- 
ative to the mass function used in the analytic calculation. The 
agreement between all three methods is excellent for masses 
below 4.2 x 10 14 M Q until our cluster catalog becomes incom- 
plete at low masses. In the most massive cluster bin, the rel- 
ative differences between the power spectra are ~ 40 - 60% 
for t ~ 2000-9000 (cf. FigH). For the lower mass bins the 
differences fluctuate between ±30%, with the pasted profiles 
generally agreeing better with the full simulation results. 

5.2. Contribution to the tSZ Power Spectrum in Redshift Bins 

In this subsection we calculate the power spectrum in red- 
shift bins and compare the results from the simulation, the 
pasted profile maps, and the analytical calculation to aid in 
understanding the differences between these approaches. In 
Figure [H] we show the power spectrum broken down into cu- 
mulative (left panel) and differential (right panel) redshift bins. 
Here we fix the mass range to M500 > 4.2 x 1O 13 M and set 
the lower redshift integration bound for the cumulative spec- 
tra to z = 0.07. We use the same definition for AQ to show 
the differences between power spectrum calculations. In con- 
trast to the mass cuts, the differences between the projected 
simulated maps and the pasted-profile maps are similar across 
all the redshift slices (cf. Fig. [8]). For £ < 5000, there is a 
~ 5 - 10% difference between the pasted profiles and the sim- 
ulations, rising to ~ 20% at £ = 10,000. These results suggests 
that the contributions from substructure and asphericity to the 
power spectrum are similar across the redshift range explored, 
with the exception of one redshift bin z ~ 0.4 which contains 
a rare merger event. The large deviations between the analytic 
and simulation/profile-paste spectra in the highest redshift bin 



are likely due to the incompleteness of the cluster catalogs at 
the lowest masses, which are preferentially more important at 
high redshift. At low redshift, we attribute the difference be- 
tween the analytic and the profile-paste power spectra to the 
shot noise in the most massive clusters. 

5.3. Contribution to the tSZ Power Spectrum within given 
Cluster Radii 

In this subsection we apply radial truncations to the full sim- 
ulated pressure maps, using clusters with M500 > 4.2 x 10 13 M Q 
at 0.07 < z < 5. The procedures for making real space radius 
cuts in maps or analytical calculations are not trivial, since any 
sharp cut in real space produces ringing in Fourier space, po- 
tentially transferring power from large to small angular scales. 
To reduce ringing and the potential to artificially increase the 
high-£ power spectrum, we use a Gaussian taper when trun- 
cating the pressure profile. We place radial tapers at r = R500, 
27?5oo, 3/?5oo, and 6R500 in the maps, adopting 6R500 as the ref- 
erence radial tapefl The form of the taper is 



T(r) = exp 



80 X flsoo 



(13) 



for r greater than the taper radius r,, and unity otherwise. In 
the bottom panel of Figure [9] we show the relative difference, 
AQ = 100 (Q,6s 500 -Q,i) /Q,6R 50 „, where Q,6s 500 is the power 
spectrum from the 6/?5oo radial cut and Q j( are power spectra 
from the other radial cuts. The trend we find is that the large 
radii of clusters are only important for the low is, for example 
the contributions to the tSZ power spectrum when only inte- 
grating out to /?5oo yields ~ 30-65% of the total power from 
I = 100- 1000, respectively. At i = 3000 only about 10% of 
the total tSZ power comes from beyond R 500. This numb er is 
consistent with previously quoted values dSun et al.112.01 ll) . We 
note that there is some small residual Fourier ringing, as the 
tapered spectra rise above the fiducial at h of many thousand. 
Nevertheless, at higher I, the cluster centers begin to be re- 
solved and become the dominant contributors to the tSZ spec- 
trum since their surface brightnesses are so much larger than 
any emission in the cluster outskirts. 

6. CONSTRAINTS OF ct 8 FROM CURRENT ACT AND SPT DATA 

Using the tSZ power spectrum and ignoring any template un- 
certainty, the constraints on 0% are competitive with other cos- 
mological measurements. After accounting for template un- 
certainty, there is no statistically significant discrepancy be- 
tween eg determined from the tSZ power and that derived 
from primary CMB anisotropics, or other t he measurements 
(IDunkley et alJ 120 lOt IShirokoff et all 1201 Oh . Here we use 
our Q.tsz templates at the fiducial parameters erg = 0.8 (and 
£lb/i = 0.03096) to define the shape of the tSZ power spec- 
trum, and content ourselves with determining only the tem- 
plate amplitude, A t sz, rel ative to that expected from the back- 
ground cosmology (e.g., iBattaglia et alJ l20Tot IDunkley et alJ 
l20"IOh . The amplit ude of A t sz is proportional to a large power 
of cr R (Atsz oc gj^lBond et alJl2002h iKomatsu & Seljakll2002l : 
iBond etalj kOOfSt iTracet all 120 111) . It follows that values of 
A t sz below unity imply that theoretical templates overestimate 

8 We avoid double counting gas particles when we project them into maps. If 
a particle lies in the overlap region between two clusters, we taper the particle 
with the larger of the two possible taper values, i.e. those particles with a 
smaller radius R/R^qq, to avoid artificially suppressing power in the overlap 
region. 




Figure 7. The tSZ power spectrum sorted into bins of cluster mass. Left: we show the cumulative tSZ power spectrum in mass bins (Q t SZ (J^soo > Moat) from the 
AGN feedback simulations, the pasted profile maps and the analytical calculation. Right: we show the differential tSZ power spectrum Q t sz (^cut.low < ^500 < 
A^cut,high) f° r me same power spectrum calculations. In the bottom of both panels we show the relative difference, AQ = 100 (Q „,„ -Q ,) /Q,j„„, where Q s , m 
is the power spectrum of the simulated maps and Qj is that from the pasted profile maps and the analytical calculation. The differences between the simulations 
and the pasted profile maps result from the absence of substructure and asphericity in the pasted profile maps, which is larger for more massive clusters. The larger 
differences found between the analytical calculation and the simulations are the resu lt of the mass catalog of the simulations having an excess of high mass clusters 
and deficit of lower mass cluster compared to the analytic mass function (cf. Fig. II U . 




Figure 8. The same as Figure [7] however for redshift slices. Left: we show the cumulative tSZ power spectrum in redshift bins Q t sz(z < Zcut) from the AGN 
feedback simulations, the pasted profile maps and the analytical calculation. Right: we show the differential tSZ power spectrum Q. t sz fcut,low < z < Zcut.high) f° r 
the same power spectrum calculations. In the bottom of both panels we show the relative difference, AQ = 100 (Q sim — Q ,) /Q „„,, where Q sjm is the power 
spectrum of the simulated maps and Qj is that from the pasted-profile maps and the analytical calculation. The agreement between the pasted profile and simulation 
spectra is excellent below £ ~ 5000 for all redshifts. On smaller scales, cluster substructure contributes similarly across all redshift bins examined. 
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Figure 9. Q t sz( r < ^cut) f° r the AGN feedback si mula tions. The thermal 
pressure distribution has been tapered as in Equation U3> at varying cluster- 
centric radii before projection. On small scales, virtually all of the power at 
t > 2000 comes from r < 2R^oo- About 80% of the tSZ power is recovered 
at I = 3000 when tapering at ftsoo, though the deviations become substantially 
larger at smaller I. These results emphasize the importance of understanding 
cluster pressure profiles well past R500 in order to do high-precision work with 
the tSZ power spectrum. 

the SZ signal, or else points to a smaller value of a% than the 
value derived from primary CMB anisotropies. 

The probability distributions of the amplitude, A t sz, and 
other cosmological parameters are determined from cur- 
rent CMB data using a modified version of CosmoMC 
dLewis & Bridle! H002), which uses Markov- Chain Monte 
Carlo techniques. We include data from WMAP 7 (lLarson et al.l 
120101) and, separately, ACT dDas et alj I20TTI) and the dusty 
star-forming galaxy-subtracted data from SPT (Shirokoff et al.1 
l2010h . We fit for 6 basic cosmological parameters (fib^ 2 > 
Odm^ 2 , «s, the primordial scalar power spectrum amplitude 
A s , the Compton depth to re-ionization r, and the angular pa- 
rameter characterizing the sound crossing distance at recom- 
bination 0) with the assumption of spatial flatness. We also 
include a white noise template for point sources Q , src with am- 
plitude A src . The pri mary difference betw een our analysis and 
the analysis by SPT ( Shirokoff et al. 2010) is that we marginal- 
ize over A src , allowing for arbitrary (positive) values, and ig- 
nore the spatial clustering component of point sources. We as- 
sume a perfect degeneracy Q.ksz oc Q jt sz for the kinetic SZ 
(kSZ) component, so we only need the relative amplitude of 
Ak,sz/A ti s z at a given frequ ency and use the kSZ amplitudes 
from Battagli a et all (1201 Oh . where the ratio o f kSZ to tSZ at 
£ = 30 00 and 150 GHz is 0.44. As mentioned in Batt aglia et al J 
d2010l) . these simulations do not fully sample the long wave- 
length tail of the velocity power spectrum a nd do not include 
any contributions from patchy re-ionization dlliev et al.l 120071 
2008). Hence this kSZ power spectrum template is a lower 
limit to the total power. 

In Figure[l0]we illustrate the 68% allowed confidence inter- 
vals for the tSZ power spectrum, given the shape of our AGN 
feedback template, our predicted tSZ-to-kSZ power spectrum 
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Figure 10. O ur 150 GHz tSZ pow er spectrum of our AGN feedback model, 
rescaled to the Keisler et al. 1 201 1) best-fit <r g value of 0.814 (red line) is con- 
traste d with the bands indicating the 6 8% range in tSZ amp litude from ACT 
IDas et alj|20Tll dark grey) and SPT jShirokoff et al. 2010, light grey). For 
comparison, we plot several other models f or the tSZ power sp ectrum, also 
shifted to the fiducial erg = 0.814. These are Sehgal et al. ( 2010) (pink dotted 
line),|Trac et al. 1 20 l j) (dark green dashed line), Komatsu & Seliak (2002) (or- 
ange dotted line), and Shaw et al. 12010) (light green dashed line). We include 
the estimated beam FWHM for ACT, SPT, and Planck. 

ratio, and the current data from ACT and SPT We scale o ur 
template using the best-fit erg value from iKeisler et al.1 (120111) of 
0.814 and scale our template (which was calculated at erg = 0.8) 
by (0.814/0.8) 8 , about 15%. We find that our template is 
within about the 68% confidence interval region for both ACT 
and SPT, after correcting for our predicted kSZ to tSZ power 
spectrum ratio of 0.44. Note that the semi-analytic and analytic 
models without substructure have lower tSZ amplitudes, which 
would result in higher values of A t sz and higher erg . 

7. DISCUSSION AND CONCLUSION 

In this work, we found a global fitting function for clus- 
ter ^herrTml_r£eMure_2rofiles using the simulations presented 
in lBattaglia et al.1 (l2010h . We find that this global fit matches 
the mean pressure profiles across mass and redshift generally 
to an accuracy of better than 10%. We have used the profile fit 
to reconstruct the thermal Sunyaev-Zel'dovich power spectrum 
using both fully analytic and semi-analytic pasted profiles onto 
cluster position in the simulations, and find we recover the tSZ 
power spectrum to ~15% at £ = 3000 (cf. Figure|5]). Other ana- 
lytic and semi-analytic models for the tSZ effect commonly as- 
sume a constant logarithmic slope of the P-p relation, T, when 
solving the equation of hydrostatic equilibrium. The assump- 
tion is not borne out in our simulations, where both the thermal 
T (which account only for the thermal pressure) and the effec- 
tive pressure T (which includes non-thermal support from bulk 
flows in clusters to the pressure) considerably increase in clus- 
ter outskirts (cf. Figure|6]l. Using both the simulations and the 
global pressure profile, we examined the contributions to the 
tSZ spectrum as functions of cluster mass, redshift, and trun- 
cation radius. We found that the contributions from substruc- 
ture and asphericity are most important for the highest mass 
clusters (M 500 > 4.2 x 10 14 M Q ), but remain significant at the 
10- 15% level across all mass bins. We find that half the power 
of the tSZ power spectrum at £ = 3000 is contributed by clusters 
with z > 0.6 and half the power originates from clusters with 
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Table 2 

Cosmological constraints on A t sz and a% from ACT and SPT using the AGN 
feedback tSZ power spectrum template 



Data 


A,sz 


°8 


ACT (Das etal. 2011) 
SPT (Shirokoff et al. 20101 


0.85 ±0.36 
0.69 ±0.29 


0.784+0;!$ 
764+0- 035 



M 5 o <2x 1O 14 M . 

We have compared our tSZ prediction to results from the 
Atacama Cosmology Telescope and the South Pole Telescope. 
We found that there is no statistically significant difference be- 
tween our model and the data, after accounting for a simplistic 
correction from the kinetic SZ effect. More complete compo- 
nent separation should be possible with better frequency cover- 
age (Millea et al. 201 lb. We note that our analysis differs from 
that in Shirok off et al.l yOlO) in that we make no prior assump- 
tion about the amplitude of the point source power spectrum, 
other than that it is non-negative. 

The pressure profile presented in this work is derived from 
the mean electron pressure in our simulations and, as such, is 
appropriate for comparison with individual clusters; we defer 
the derivation of a mean profile designed to include the effects 
of substructure and asphericity in the power spectrum to a fu- 
ture work (Battaglia et al., in prep). This profile will not be 
expected to match individual cluster observations, but we hope 
will allow analytic calculations of the tSZ power spectrum to 
an accuracy of significantly better than 10%. With future data 
sets, such as those expected from Planck, ACTpol, and SPT- 
pol, it may be possible to constrain not just the amplitude but 
the shape of the tSZ spectrum. In this case, analytic calcula- 
tions may be usable to constrain not just cosmology but the 
important astrophysical processes in clusters with the tSZ ef- 
fect. Doing so through the power spectrum has the advantage 
that it is sensitive to lower mass and higher redshift clusters 
as well as cluster outskirts in ways that are complementary to 
other data sets. 
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Figure 11. We compare the mass function, An/ AM, fo r the cluster catalog 
from AGN feedback simulations to the mass function from Tinker et al. 12008). 
The differences at high masses indicates that in the 10 independent simula- 
tions we happen to have more high-mass clusters than is expected on average 
(though with only 6 with M500 > 7.1 X 10 14 Mq this is consistent with shot 
noise). At low masses, our catalog is incomplete due to our FOF halo finding 
(see text). 

APPENDIX 

A. COMPARING THE CLUSTER MASS CATALOG TO THE MASS 

FUNCTION 

In this appendix we compare the mass function from our 
simulations with that of iTinker et al.l (120081) . Our cluster mass 
catalogs were made with spherical overdensity mass with re- 
spect to the critical density and the mass function is with re- 
spect to the mean matter density. So, we converted the M200 
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from the simulations to M2001B assuming the mass profile is 
dominated b y dark matter and u se the concentration-mass re- 
lations from iDuffv et al.l (120081) . We show in Figure [TT1 that 
there is a clear deficit of low mass clusters due to the chosen 
linking length of 0.2 in our FOF finder. At this length, is is 
well known that neigh boring clusters are sometimes artificially 
merged together (e.g.. iDavis et al.lfl985l) . We also instituted 
a firm lower limit mass cutoff in the initial FOF catalogs of 
M F0F > 1.4 x 10 I3 M Q , and so our mass function is also ex- 
pected to be incomplete near that mass. 

There is a clear excess of high-mass clusters in our simu- 
lations, but it is consistent with shot noise (we only have 6 
clusters with M500 > 7.1 x 10 I4 M Q ). We now estimate the ex- 
cess power in our full simulation power spectrum due to this 
upwards fluctuation in the highest mass bin. Where the cluster 
catalogs are complete, we expect that over an enormous num- 
ber of simulations, the paste profile and analytic calculation 
of the tSZ power spectrum would converge, and indeed see 
the agreement is excellent between the two in the right panel 
of Figure [7] for all but the lowest (due to catalog incomplete- 
ness) and highest (due to shot noise) mass bins. We therefore 
adopt the ratio of the the pasted profile spectrum to the analytic 
spectrum as a quantitative estimate of the over-representation 
of high mass clusters in our finite number of realizations. At 
t = 3000 this ratio is 2.0 for clusters with M 50 o > 7. 1 x 1O 14 M , 
though as can be seen in Figure |7]the specific value is insensi- 
tive to the reference £. Since the high-mass contribution to the 
tSZ spectrum from the full simulation projections is 0.67/iK 2 
at I = 3000, in the limit of an infinite number of simulations, 
we would expect the average contribution from clusters with 
M 500 > 7.1 x 10 14 M Q to be 0.34/xK 2 lower. The total power 
spectrum at I = 3000 is 5.78 /iK 2 , so this shot noise correction 
amounts to just less than a 6% shift in the total power spectrum. 



